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Abstract 

The CAT Imaging Telescope has observed the BL Lac object Markarian 501 between 
March and August 1997. We report here on the variability over this time including several 
large flares. We present also preliminary spectra for all these data, for the low emission 
state, and for the largest flare. 



Cj " 1. Introduction and Data Sample 

o' 

Several groups have reported strong emission from Mrk501 at this workshop. The CAT Imaging 
■ Telescope, operating since October 1996 and (described in Rivoal et al. 1997), is located in the 
French Pyrenees (42,5° N, 2° E) at an altitude of 1650 m. It has a mirror area of 17 m 2 with a 
fine grained imaging camera (546 pixels, each 0.12° for these data) and a threshold energy of 
OS | 220 GeV. Between March and August 1997 a total of 80 hours of data were taken on Mrk501 
with 25 hours of control region (OFF) data. As the performance of the telescope at large zenith 
angles is still under investigation, the data used for the variability study were limited to those 



at less than 25° from Zenith, and for the spectral analysis at less than 10° (16.58 hours). 
2. Analysis method 
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A standard moment-based ( "SuperCuts" ) analysis has been used for the study of Mrk501's 
light curve (section 3). The cuts used are lower than for other detectors, taking advantage of 
the CAT telescope's very high resolution imaging camera, as follows: 
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I. 

A new method has been developed for the CAT telescope (Lebohec 1996 & Lebohec et 
al. 1995) which provides directly the energy of the incident 7-ray; this has been used in the 
spectral analysis. The method uses a maximum-likelihood fit to the images of an analytical 
model for 7-shower images. The image is given by the model as a function of primary energy, 
impact parameter (R), angular source position in the field of view and azimuthal position of the 
image about the source. Minimization of the x 2 fit provides an estimate of these parameters. 
In addition, the minimized \ 2 value gives a criterion for selection of good 7-ray candidates. 
The orientation angle a in this method is defined (between and 180°) as the angle at the 
image centre between the expected source position and reconstructed source position. Figure 1 
shows the clear signal in the alpha plot of the data from Mrk501 at less than 25° from Zenith. 
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Fig. 1: On-source and Off-source a- distributions for March to August observations of Mrk501. 

The origin of each gamma shower on the celestial sphere can be reconstructed event by 
event to 0.16°. The energy resolution obtained with this method is 20%, independent of energy 
from 200 GeV to 10 TeV. 

Three cuts are used for the spectral analysis (derived from simulations to optimize the 
number of 7-ray candidates): 

0.2 < P X 2 

40.0 < R (metres) 
a < 9 (deg.) 



3. Source Variability and Light Curve 

The nightly flux levels (expressed in gamma-rays per minute uncorrected for cut efficiencies) 
for Mrk501 between March 7 and August 8 1997 are shown in figure 2. During this period 
the average flux was higher than that of the Crab Nebula, with the highest flare over 6 times 
greater. Variability on a day time-scale is evident in this figure; shorter time variations are still 
under study. 

Figure 3 shows the remarkable agreement between the data from three observatories 
(CAT, Whipple, and HEGRA) for the April period. The data from each experiment have been 
rescaled as they have different thresholds and observatory altitudes. This demonstrates also the 
utility/necessity of having several observatories to follow continuously such variable sources. 
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Fig. 2: Nightly gamma-ray rate for Mrk501 from March to August, 1997. Coverage of the source 
was reduced in May due to technical difficulties, and in August due to weather conditions. 



4. Differential Energy Spectra 

The acceptance area and the cut efficiency were derived from detailed Monte Carlo simulations 
as a function of energy and zenith angle. These, combined with the energies from the maximum- 
likelihood fit, allow differential energy spectra to be obtained. Spectra from Mrk501 are shown 
in figure 4 (for data less than 10° from Zenith): for all data, for the data for which the source 
was in a low state, and for the largest flare on the night of April 15-16th. A maximum-likelihood 
fit of a power law to the spectrum of these data gives for all data: 
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The corresponding flux and index for the low state data are 2.50 ± 0.18 and —2.59 ± 0.09, and 
for the largest flare 18.94 ± 0.64 and -2.33 ± 0.05. 

The errors in these formula are statistical only as the systematics studies are underway. 
Within statistical errors the spectrum index seems compatible with a constant hardness of 
the source. A cross-check of these values by the standard SuperCuts analysis yields closely 
compatible results both for flux values and spectral indices. The spectra appear to be harder 
than that of the Crab Nebula (given in Goret et al., 1997). In this dataset alone, there is a 
4.7<t signal at energies above 7 TeV; if the cut on P x 2 is relaxed the signal increases to 6.3cr 
above 7 TeV, with 4.8a above 10 TeV. Further studies using the data taken at high Zenith 
angle should extend the spectral coverage to higher energies. 
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Fig. 3: April 1997 nightly gamma-ray rate for Mrk501 as seen by three groups : CAT, Whipple 
(Fegan, 1997), and HEGRA system (Lorenz, 1997) . 



5. Conclusions 

The source Markarian 501 entered a very active phase in 1997, with several large flares during 
which it became the brightest object in the VHE 7-ray sky. This unusual behaviour has 
been detected by the CAT imaging telescope and several other groups. The intensity of the 
emission provides data which are exteremely rich in 7-rays, useful both for testing of the detector 
capabilities and investigation of the highest-energy part of the spectrum. At the time of writing 
the source has been seen to be still active in data taken at high Zenith angle. 
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Fig. 4: Differential spectra for Mrk501 for all data, low-state emission, and the highest Bare. 



